EXISTENCE, UNIQUENESS AND SMOOTHNESS OF SOLUTION FOR 
3D NAVIER-STOKES EQUATIONS WITH ANY SMOOTH INITIAL VELOCITY. 
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Abstract. Different authors had received a iot of results regarding the Euler and Navier-Stokes equations. Existence and 
smoothness of solution for the Navier-Stokes equations in two dimensions have been known for a long time. Leray [3] showed 
that the Navier-Stokes equations in three dimensional space have a weak solution. Scheffer [4], [5] and Shnirelman [6] obtained 
weak solution of the Euler equations with compact support in spacetime. Caflfarelli, Kohn and Nirenberg [J improved Schcffcr's 
results, and F.-H. Lin [8] simplified the proof of the results of J. Leray. Many problems and conjectures about behavior of weak 
solutions of the Euler and Navier-Stokes equations are described in the books of Ladyzhenskaya [I], Bertozzi and Majda Q5], 
O | Tcmam |10l . Constantin 1 1 1 1 or Lemari-Rieusset 1121 . 

Solutions of the Navier-Stokes and Euler equations with initial conditions (Cauchy problem) for 2D and 3D cases were 
i obtained in the converging series form by analytical iterative method using Fourier and Laplace transforms in paper 1 1 41 . These 

solutions were received as infinitely differentiable functions. That allowed us to analyze essential aspects of the problem on a 
much deeper level and with more details. For several combinations of problem parameters numerical results were obtained and 
presented as graphs |14| . 1151 . 

This paper describes detailed proof of existence and uniqueness of the solution of the Cauchy problem for the 3D Navier- 
Stokes equations with any smooth initial velocity. This solution satisfies the conditions required in I13| for the problem of 
Navier-Stokes equations. When viscosity tends to zero this proof is correct for the Euler equations also. 

1. Introduction. 

The solution of the Cauchy problem for the 3D Navier-Stokes equations described in this paper is based 
on the form of differential equations in the statement of the problem and also conditions for the initial 
velocity and applied force [13]. It grows from classic definition of function and classic methods of analysis. 

In the problem statement for Navier-Stokes equations an initial velocity is infinitely differentiable function 
decreasing rapidly to zero in infinity. Applied force is identically zero. 
Solution of the problem will be presented by the following stages: 

First stage (sections 2, 3) - we have moved non-linear parts of equations to the right sides. Then we 
have solved the system of linear partial differential equations with constant coefficients. We have obtained 
the solution of this system using Fourier transforms for the space coordinates and Laplace transform for time. 
From theorems about application of Fourier and Laplace transforms for system of linear partial differential 
equations with constant coefficients we see that in case if initial velocity and applied force are smooth enough 
functions decreasing in infinity, then the solution of such system is also a smooth function. (Corresponding 
theorems are presented in S. Bochner [IB], V.P. Palamodov [T7|, G.E. Shilov [TH], L. Hormander [TH], S. 
Mizohata [20], J.F. Treves [H]). Result of this stage is the integral equation for vector-function of velocity. 

Second stage (sections 4, 5) - we have introduced perfect spaces of functions and vector-functions (I.M. 
Gel'fand, G.E. Chilov [22]): in which we have looked for the solution of the problem. We have demonstrated 
equivalence of the solution of the Cauchy problem in forms of differential and integral equations. 

Third stage (section 6) - we have divided all parts of integral equation by an appropriate constant 
depending on value of initial fluid velocity, and received equivalent integral equation. We have also cor- 
respondingly replaced integration variables in integral operators. This newly received equivalent integral 
equation allowed us to analyze the Cauchy problem for the 3D Navier-Stokes equations for any value of 
initial fluid velocity. According to a priori estimation of the solution of the Cauchy problem for the 3D 
Navier-Stokes equations [T], [2] the described constant is proportional to max of the norms of the initial 
velocity in spaces C 2 and L 2 . 

Fourth stage (section 6) - We have used the newly received equivalent integral equation to prove exis- 
tence and uniqueness of the solution of the Cauchy problem in all time range [0,oo) based on the Caccioppoli- 
Banach fixed point principle (L.V. Kantorovich, G.P. Akilov [23], V.A. Trenogin [23, W. Rudin [25], W.A. 
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Kirk and B. Sims [55] , A. Granas and J. Dugundji ^7\, J.M. Ayerbe Toledano, T. Dominguez Benavides, G. 
Lopez Acedo |28j ). For this purpose three following theorems were proven in this paper: 
Theorem 1: Integral operator of the problem is a contraction operator; 

Theorem 2: Existence and uniqueness of the solution of the problem is valid for any t G [0,oo); 
Theorem 3: Solution of the problem is depending on t continuously. 

Fifth stage (section 6) - By using a priori estimation of the solution of the Cauchy problem for the 3D 
Navier-Stokes equations [T|, P] we have shown that the energy of the whole process has a finite value for 
any t G [0,oo). 

2. The mathematical setup. 

The Navier-Stokes equations describe the motion of a fluid in R N (N = 3). We look for a viscous 
incompressible fluid filling all of R N here. The Navier-Stokes equations are then given by 

(2-1) ^ + = ^ - ^ + AOM) (x£R N , t > 0, 1 < k < N) 



(2.2) dwu = V ^ = (xeR N ,t>0) 

, ox„ 

n—l 



with initial conditions 



(2.3) u(x,0) = if(x) (x G i? JV ) 

Here u(x,t) — (uk(x,t)) G R N , (1 < k < N) — is an unknown velocity vector (N — 3); p(x,t) — 
is an unknown pressure; v?(x) is a given, C°° divergence-free vector field; fk{x,t) are components of a 
given, externally applied force f(x,t); v is a positive coefficient of the viscosity (if v = then (|2.1[) - (|2.3p 
are the Euler equations); and A = J^Li is the Laplacian in the space variables. Equation (|2.ip is 
Newton's law for a fluid element subject. Equation (|2.2|) says that the fluid is incompressible. For physically 
reasonable solutions, we accept 

(2.4) u k (x,t) , ^ as I x I -> oo (1 < k < N, 1 < n < N) 
Hence, we will restrict attention to initial conditions m° and force / that satisfy 

(2.5) | d"y°(x) | < C aK (l+ | x \)~ K on R N for any a and K. 
and 

(2.6) | d"d?f{x, t) | < C a0K (l+\x\ +t)- K on R N x [0, oo) for any a, P and K. 

To start the process of solution let us add (— u nTfT~ ) ^° both sides of the equations (|2.1[) . Then 

we have: 
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(2.7) ^ = v^u k - $L + f k (x,t)- f>f^ (x€R N , t>0, l<k<N) 

n—1 n 

N r. 

(2.8) divu = V = (xeR N ,t>0) 

(2.9) u{x,0) = vP{x) (x G R N ) 



(2.10) u k (x,t)^0 , ^ as I x I -> oo (1 < k < N, 1 < n < N) 



(2.11) | 9£u°(a;) | < C7 Qif (1+ | x \y K on R N for any a and K. 

(2.12) | d%d?f(x,t) | < C Q(3K (1+ | x | +t)~ K on R N x [0, oo) for any a,/3 and K. 
Let us denote 



(2.13) / fc (a:,f) = f k (x,t) - ^«„^ (1 < k < N) 

n—1 



or we can present it in the vector form: 
(2.14) f(x,t) = f(x,t) - (u ■ V)u 



3. Solution of the system of equations ()2.7p - (|2.14|) . 

Let us assume that all operations below are valid. The validity of these operations will be proved in the 
next sections. Taking into account our substitution f|2 . 13[) we see that equations (|2.7p — (|2.9|) are in fact 
system of linear partial differential equations with constant coefficients. 

Solution of this system will be presented by the following steps: 

On the first step of our analysis we use Fourier transform (|A1[) to solve equations (|2.7p — (|2.14p . We 
have got: 



#fc (71j72>73i*) = F[u k (x 1 ,x 2 ,x 3 ,t)} 



2tt i +\ v\ d 2 u k (x 1 ,x 2 ,x 3 ,t) ^ 
-7 s C / fc(7i)72,73,*) = F[ J [use p.lfjp] 
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#£(71,72,73) = F[ul(xi,x 2 ,x 3 )] 
^(71,72,73,0 = F\p(x 1 ,x 2 ,x 3i t)] 
F k (ji,l2,l3,t) = F[f k (x 1 ,x 2 ,x 3 ,t)} 
k,s = 1,2,3 



and then: 



(3-1) ^^^P' 73 '^ = +72 + 7i)#i (7i, 72, 73,*) + i7i-P(7i,72,73,*) + Fi (71, 72, 73,0 



^#2(71, 72, 73,0 
dt 



(3-2) "r - " — — = -^(7i +72 +73)^2(71, 72, 73,0 + «72-P(7i,72,73,*) + ^2(71, 72, 73, 



(3-3) rfC/3 ^ 7l ^ 2 ' 73 ' f ^ = -^(7i +72 +73)^3(71, 72, 73,0 + «73-P(7i,72,73,0 + ^3(71, 72, 73,0 



(3.4) 7i #1(71, 72,73,0 + 72 #2(71,72,73,0 + 73 #3(71,72,73,0 = 

(3-5) #i(7i,72,73,0) = #1(71,72.73) 

(3-6) #2(71,72,73,0) = #2(71,72,73) 

(3-7) #3(71,72,73,0) = #3(71,72,73) 

Hence, we have received a system of linear ordinary differential equations with constant coefficients 
(|3.1[) — (|3.7p according to Fourier transforms. At the same time the initial conditions are set only for 
Fourier transforms of velocity components U\ (71 , 72 , 73 , , #2 (71 , 72 , 73 , , #3 (71 , 72 , 73 , ■ Because of that 
we can eliminate Fourier tranform for pressure P(7i, 72, 73, from equations p.lj) — (|3.3|) on the second 
step of solution. From here assuming that 71 7^ 0, 72 7^ 0, 73 7^ 0, we eliminate P("fi, 72, 73, from 
equations (|3.ip — (|3.3p and find: 

^t[#2(7i,72,73,0 - —#i(7i,72,73,0] = ~f(7i +72 +7f)[ #2(71,72,73,0 ~ 
at 71 

(3-8) -—#i(7i,72,73,0] + [-^(71,72,73,0 ~ —-Pi(7i,72,73,0] 

7i 7i 
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-7i [ U 3 (ji, 72,73, *) -i7i(7i,7a,73)*)] = - v ill +ll +lt)[ ^3(71, 72, 73,*) - 



(3.9) -^i(7i,72,73,*)] + [ -£3(71,72,73,*) --Pi(7i,72,73,*)] 

7i 7i 

(3.10) 7i^i(7i,72,73,*) + 72^2(71,72,73,*) + 73^3(71,72,73,*) = 

(3.11) f^i(7i,72,73,0) = ^"(71,72,73) 

(3.12) C/ 2 (7i,72,73,0) = t/ 2 °(7i,72,73) 
(3-13) C^3(7i,72,73,0) = 1/3(71,72,73) 



On the third step we use Laplace transform (|A2j) . (|A3j) for a system of linear ordinary differential 
equations with constant coefficients (|3.8p — (|3.10[) and have as a result the system of linear algebraic 
equations with constant coefficients: 

U® (71, 72,73, v) = ^[^(71,72,73,*)] k=l,2,3 
^f(7i,72,73,77) - L[F fe ( 7l ,72,73,*)] k=l,2,3 

v[ Uf (71, 72,73,?7) - —Uf (71,72,73,7?)] - [^2(71,72,73,0) - — £/i(7i,72,73,0)] = 
7i 7i 

+ 72 +73)[ U f (71,72,73,77) ~ —^(71,72,73,77)] + 

7i 

(3-14) +[^(71,72,73,77) - --Ff (71,72,73,77)] 

7i 

7? [ Uf (71, 72, 73, 7?) - —Uf (71,72,73,77)] - [^3(71,72,73,0) - — I7l(7l,72.73,0)] = 
7i 7i 

+72 ^(71,72,73,77) - — Uf (71,72,73,7?)] + 

7i 

(3-15) + [ Ff (71, 72, 73, 7?) - -Ff (71,72,73,77)] 

7i 

(3.16) 71^(71,72,73,77) + 72^(71,72,73,77) + 73 £/f (71,72,73,77) = 
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(3.17) 



C^i (7i)72, 73.0) = ^1(71,72,73) 



(3.18) 



^2(71,72,73,0) = ^2(71,72,73) 



(3-19) ^3(71,72,73,0) = ^"(71,72,73) 

Let us rewrite system of equations (|3 . 14[) — (|3.16[) in the following form: 



V + ^(7i 2 + 7 2 2 + 7 3 2 )] —U? (71,72,73,??) - [v + "(Ti + 7 2 2 + 7 3 2 )] ^2®(7i, 72, 7s, rj) 
7i 



(3.20) [ — Ff(7i, 72,73, V) ~ F® (71)72,73,??)] +[ — fi(7i, 72, 73, 0) - ^2(71, 72, 73,0)] 
7i 7i 



V + + 7 2 2 + 7 3 2 )] jU? (71,72,73,^) - [v + + 7 2 2 + 7 3 2 )] ^(71,72, 73, rj) 



(3.21) [ —Ff (71, 72,73, rj) - Ff (71,72,73,^)] + [ — C^i(7i,72,73,0) - t/3(7i,72,73,0)] 
7i 7i 



73, 



(3.22) 71 Uf (71,72,73, »7) + 72 Uf{ lx , 72,73, r?) + 73 ^f(7i, 72,73,7?) = 

Determinant of this system is 



A 



V + K7i 2 +72 2 +7 3 2 )] £ 
r? + K7i 2 +7 2 2 +7 2 )] ^ 
7i 



-[ V + K7i 2 +72 2 +7f; 



02 





v + +72 +73) 
73 



(3.23) 



7/ + K7i + 7 2 2 + 7 3 2 )] (7i 2 + 7 2 2 + 7 3 2 ) 
7i 



^ 



And consequently the system of equations (|3.14|) — (|3.16| and/or (|3.20[1 — (|3.22| has a unique 
solution. Taking into account formulas (|3.17| — (I3.19| we can write this solution in the following form: 



uf (71,72,73,7?) 



(3.24) 



[(72 + 7 2 ) F f (7i , 72 , 73 , rj) ~ 7i72 Ff (71 , 72 , 73 , y) - Tj/yjsjf (71 , 72 , 73 , v)] 
(7? + 7 2 2 + + " (7i + ll + 7 3 2 )] 

^1(71,72,73) 



[V + f(7i + ll + 7 3 2 )] 
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^(71,72,73,77) 



(3.25) 



[(73 + 7i)-Pf (7i, 72, 73, ??) - 7273-Ff (7i, 72, 73, V) ~ 727i Ff (71, 72, 73, ??)] 
(7? + 7| + 7s) [77 + K7i 2 + 7| + 7 3 2 )] 

^2(71,72,73) 



[V + v(ci\ + ll + 7 3 2 )] 



^(71,72,73,77) = 



(3.26) 



[(li + li ) Ff (7i , 72 , 73 , V) ~ 737i Ff (7i , 72 , 73 , 77) ~ 7372 Ff (71 , 72 , 73 , y)] 
(7i + 7l + 7 3 2 ) N + K7 2 + 7 2 2 + 7 3 2 )] 

^3(71,72,73) 



b? + K7i + 7l + 73)] 



Then we use the convolution theorem with the convolution formula (|A.4[) and integral (|A.5[) for (|3.24p 
(|3.26p and obtain: 



^i(7i,72,73,*) 



_ i , (7 2 +7 2 +7 2 )(t _ T) [(71 + 71)^1(71,72,73,7-) -7172^2(71, 72, 73, r) -7i73-F3(7i,72,73,t)] ^ 

(7? + 7l + 7 3 2 ) 

(3.27) + e -"(T?-Hrf+73 a )t C/0( 7l )72j73 ) 



^2(71,72,73,*) 



/o 
(3.28) 



e -^(7?+7, 2 +7^)(t-r) [(73 +7i 2 )A(7i,72,73,t) - 7 2 73-F 3 (71 , 72 , 73 , r) - 7271-^1(71, 72, 73, r)] 

(7i 2 + 72 2 + 7 3 2 ) 

+ e -W+^+^ 2 )'c7 2 o (7l)72!73) 



^3(71,72,73,*) 



/() 
(3.29) 



P ->A^+ilW,){t-r) [(7i+ 72)^3(71, 72,73, t) -7371-^1(71, 72, 73, t) - 7372-^2(71, 72, 73, r)] 

(7! 2 +7 2 2 + 7 3 2 ) 

+ e -^ 2+ ^ + ^ 2)t t/ 3 (7i,72,73) 



Using the Fourier inversion formula (|^4. 1|) we receive: 



u 1 (x 1 ,x 2 ,x? n t) 



(2tt) 3 / 2 



OO pOO /•DO 



00 ./ — 00 «/ —00 



e -<^+~fZ+~fi)(t-r) [(7 2 2 + 7 3 2 )A(7i,72,73,r)] ^ 

(7i 2 +72 2 + 7 3 2 ) 
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/"* e - i ,( 7 ? +7 ?+7, 2 )(t-T) [7172^2(71, 72, 73, t) +7173^3(71: 72,73, r)} ^ 
J a (7i 2 +7 2 2 +7|) 



+ e -^ + ^ + ^)* [/ o (71j72j7 3 ) 



1 

8^ 



1 

8^ 



1 

8^ 



OO /> OO pOO 



OO </ — OO J — OO 



00 poo />oo 



(7I + 7I) 

(7i 2 +7 2 2 +7 3 2 ) 



g _ 1/ ( 7 3 +7 2 +7 2 )(t _ r) 



OO /"OO /»oo 



— 00 •/ — 00 «/ — OO 



i(5i7i +£272+^373) 



— 00 J — 00 



7172 



— 00 (7i 2 + 7 2 2 + 73) 



■fi{x\, X2, X3,t) dx.1dx2dx.3dr 



g -f(7i +75+7 3 )(*~ T ) 



g -i( xl 7 1 +x 272 +X37 3 ) d 7l rf 72 ^ 73 



OO />00 /"OO 



— 00 J —00 J — OO 



71+^272+^373) 



00 /> 00 /< 00 



00 </ —00 </ — OO 



7173 



•fzixi, X2, X3,t) dxidx 2 dx 3 dT 



e _^( 7l ^ +7 ^ +7 ^)(t_ T ) 



e -i(x 17 i+x 2 72+X373) d 7l d 72 d 73 



(7i 2 +7 2 2 +7l) 

•/ 3 (ii, i 2 , £3, t) dxidx 2 dx 3 dT 



00 />oo /"OO 



— 00 J —00 J —00 



i(xi 71 +£272+^373) 



e -*(^lTl+^272+X373) d7ld72d73 + 



+ 



8tt 3 



OO /"OO /"OO 



— 00 </ — 00 ./ — 00 



g -"(7f+Tj+7|)< 



00 />oo />oo 



— 00 ./ — 00 ./ — 00 



j(xi7i+£272+X37 3 ) 



iti(xi, i 2 , i 3 ) dx\dx2dxz 



g - t (x 1 7 1 +X 27 2+X373) rf 7l rf 72 d 73 



(3.30) 



ftiC/i) + Sia(/ 2 ) + Si3(/ 3 ) + B(«S), (1231) 



w 2 (xi,a; 2 ,x 3 ,t) = 



(2tt) 3 / 2 



00 />oo />oo 



OO ^ — 00 ^—00 



1) 



-vh?+^+^)(t-r) [hi +7l)^2(7l, 72,73, t)} ^ _ 

(7i 2 + 7 2 2 + 7 3 2 ) 



e - 1 ^( 7l 2 +7, 2 +7, 2 )(t-T) [7273-^3(71, 72, 73, t) + 7271.F1 (71, 72,73, r)] 
a (7i 2 +7 2 2 + 7 3 2 ) 



+ e -^ + ^ + ^)*c/ 2 ( 7l ,7 2j 7 3 ) 



g -i(x l7 i+x 2 7 2 +X373) d 7l d 72 d 73 



1 

8^3 



OO /"OO />00 



— 00 </ — OO 



7271 



g -" y (7i+72+7 3 )( t — r ) 



-00 (7i + 7 2 2 + 7 3 2 ) 

■fi (x 1 , x 2 , £3, t) dxidx 2 dx 3 dT 



00 poo poo 



— 00 J —00 J —00 



i(xi7i+x 2 72+x 3 7 3 ) _ 
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8tt 3 



— oo J — oo 



-oo (7? + ll + ll) 

■f2{x\,X2,X3,r) dxidx2dx 3 dT 



oo />oo f oo 



i(5i7i+X272+i373) 



1 

8^ 



— oo J — oo 



7273 



g i(xi7l +£272+2373) 



6 -K7r+72+7 3 )(*-^) 



-oo (7? + ll + ll) 

•f3(xi, X2, X3, t) dx\dX2dX3,dT 



— oo J — oo <J — oo 



oo />oo />oo 



— oo — oo — oo 



g -i (a:i7l+a;272+ x373) d 7l d 72 d 73 + 



1 

8^ 



OO />0O />OC 



oo J -co J — oo 



e _„( 71 - +7 2 +7 - )t 



oo />oo /"OO 



— oo J — oo J — oo 



i(£i7i+2 2 72+£ 3 73) 



u%{£\, £2, X3) dxxdx^dx^ 



e ~ l{Xi ' <i+X2 ' <2+X3 ' t3 U~f 1 d l2 d> l3 = 



(3.31) 



5 2 i(/i) + S 2 2{~h) + S 23 (h) + B(u° 2 ), (AM 



u 3 (xi,x 2 ,x 3 ,t) 



(2tt) 3 / 2 



OO /-OO {-CO 



00 J —00 J —00 



M^+jj+JsKt-r) [(7i +72)^3(71. 72,73, r)] ^ 
(7i 2 + 7 2 2 + 7 3 2 ) 



-^(7! 2 +7j+7g)(t-T) [7371-^1(71, 72, 73, T) +7372-^2(71, 72,73, T)] 



(It 



+ e -(7 1 2 +72 2 +7 3 2 )t £/ 3 °( 7l , 72 , 73 ) 



(7i 2 +7 2 2 + 7 3 2 ) 

g - J :(x l7l + 2; 272+^373) rf 7l d 72 d 73 



1 

8^3 



1 

8^ 



OO /*0O /'OO 



7371 



00 ./ —00 J —00 



OO /'OO /*00 



(7i 2 + 7 2 2 + 7 3 2 ) 



g _„( 7 2 +7 2 +7 2)( t _ T ) 



00 /*oo /*oo 



00 «/ —00 J —00 



Q 1(^171+^272+^373) 



#2, ^3? t) dx\dx2dx^dr 



00 «/ —00 «/ —00 



7372 



(7i 2 +7| + 7 3 2 ) 



j_ r /-°° r _(7i 2 +7 2 2 ) 

8tt 3 



00 «/ —00 «/ —00 



(7i 2 + 7 2 2 + 7 3 2 ) 



/OO r OO /*00 
Q i(^l7l+^272+^373) 
-00 */ —00 */ —00 

/•OO />00 />00 

g -^(7?+72+7 3 )( t ~ T ) / / / e i (*171+»a72+*373) 



• /2 (i 1 , x 2 , x 3 , t ) dx 1 dx 2 dx 3 dr 



00 J — 00 — 00 



■f3(xi,X2,X3,T) dxidx 2 dx 3 dT 



g -i(*l7l+<l*7a+*3T*) dlldl2dl3 + 



1 

8^3 



OO /*oo /«OC 



00 ./ —00 «/ — OO 



£ -^(7r+7 2 +7 3 )* 



OO pOO /*oo 



CO «/ — OO ./ —CO 



q i(xiji +^272+2:373) . 
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• u° 3 (x 1 ,x 2 ,x 3 )dx 1 dx 2 dx 3 e - li - Xl ^ +X2 ^ +x ^ d lx d l2 di 3 = 

(3.32) = S-ai(A) + S 32 (h) + S 33 (f 3 ) + B(u° 3 ), (031 

HereSuQ, S i2 Q, S 13 (), S 21 (), S 22 Q, S 23 (), S 31 {), S 32 Q, S 33 {), BQ are the integral operators. 

S 12 Q = <M) 
SisQ = S31O 
S 23 Q = S 32 Q 

It follows from (|3.30[) — (|3.32|) that the vector u can be represented as: 

(3.33) u = S- f + B-^ = S-f-S-(u-V)u + B-u°, 

where / is determined by formula (|2.14p . 

Here S and B arc the matrix integral operators: 



Sn 


S12 


<Sl3 




5*22 


S23 


5*31 


5*32 


S33 



B 











B 











B 



4. Spaces S, TS. [22], [29] 



JV 



Let us consider space S of all infinitely diffcrcntiable functions <p(x) defined in N-dimcnsional space R 
(N = 3), such that when | x | — > 00 these functions tend to 0, as well as their derivatives of any order 
more rapidly than any power of t^t. 

To define topology in the space S let us introduce countable system of norms 

(4.1) y\\ p = su Px \x k D«y(x)\ (p = 0,1,2,...) 



where 



JJt P(X)\ — \x 1 ...x N „ qi — jr-m— 



k = (ki,...,ktf), q= (qi,...,q N ), x k = x^ 1 . . . 
D " = DJ'.'.'.Z^ ■ (fci,---, 9^ = 0,1,2,...) 
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Space S is a perfect space (complete countably normed space, in which the bounded sets are compact). 
Space TZ> of vector-functions <p is a direct sum of N perfect spaces S (N = 3) [24], ED] : 

TS = S®S®S. 

To define topology in the space let us introduce countable system of norms 

N N 

(4.2) \\0\\ p = YsWViWp =J2 SU P* \x k D q Vi(x)\ (p = 0, 1, 2, ...), (N = 3) 

i= i i= i 1*1, kl <p 

The Fourier transform maps the space S onto the whole space S and maps the space TS~ onto the whole 
space T§ [IS], [22]. 

5. Equivalence of Cauchy problem in differential form (12 . 1|) - (|2.3p and in the form of an 
integral equation. 

Let us denote solution of the problem (|2.1[) - (|2.3|) as {it(xi, X2, X3, i), p(xi, %2, X3, t)}, in other words 
let us consider infinitely differcntiablc by t € [0,oo) vector-function it(xi, X2, X3, t) e TS and infinitely 
differcntiablc function p(xi,X2,X3, t) G S, that turn equations (|2.1[) , (|2.2j) into identities. Vector-function 
u(xi, x 2 , X3, i) also satisfies the initial condition (|2.3[) (u°(xi, X2, X3) 6 TS): 

(5.1) u(xi,x 2 ,x 3 ,<)|t =0 = u°(xi,x 2 ,x 3 ) 

Let us put {u(xi, X2, X3, t), p(xi,X2,X3, t)} into equations (|2.ip , (|2.2[) and apply Fourier and Laplace 
transforms to the result identities considering initial condition (|2.3|) . After all required operations (as in 
sections 2 and 3) we receive that vector- function u(xi, X2, X3, t) satisfies integral equation: 

(5.2) u = § ■ f -B-(u- V)u + B ■ m° = 5 V • u 

Then vector-function grad p £ TS is defined by equations (|2.ip where vector-function u is received from 
equation (|5.2[) . 

Here / G T^, it G T^ 1 and S, B, S are matrix integral operators. 

Vector-functions S ■ /, B • it , S ■ (u ■ V)u are also belong T5 since Fourier transform maps perfect 
space onto TS. 

Going from the other side, let us assume that u(x\, x 2 , X3, t) G TS is continuous in t G [0,co) solution 
of integral equation (|5.2p . Integral-operators Sij ■ (u • V)u are continuous in t G [0,oo) [see p.30p — p.32|) ]. 
From here we receive that according to (|5.2[) 

u(xi,x 2 ,x 3 ,0) = 72°(xi,x 2 ,x 3 ) 

and also that u(x\, x 2 , X3, t) is differcntiablc by t G [0,co). As described before, the Fourier transform 
maps perfect space TS on itself. Hence, {w(xi, X2, X3, t) and p(xi, X2, X3, t)} is the solution of the Cauchy 
problem (|2.1[) - (|2.3p . From here we see that solving the Cauchy problem (|2.1[) - (|2.3p is equivalent to finding 
continuous in t G [0,oo) solution of integral equation (|5.2[) . 
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6. The Caccioppoli-Banach fixed point principle. [23], [23], [25], [25], [27], [25] 

From the problem statement [13] we have / = 0. Let us rewrite integral equation (|5.2p with this 
condition 

(6.1) u = -5- (u- V)u + 

Let us divide all parts of the integral equation (16. ip by some constant V, that we will define appropriately 
below. Then we receive modified integral equation equivalent to equation (|6.1|) : 



(6.2) u v = —S v ■ [Uy ■ V v )u v + By ■ u v — S v -m v 
here 

— -o 

(6.3) u v = -; 4 = — ; V v = V ■ V; x kv = 7fcv = V ■ 7fe ; ^ v = ^ ; (1 < fe < N) 

The selection of constant V is based on a priori estimation of the solution of the Cauchy problem for 
the Navier- Stokes equations p], [2j: 

(6-4) \\u\\ L2 < |K|| i2 

where 

. V2 

|l 2 = f f R3 \u\ 2 dx 

The space of vector - functions T$ is divided into two subsets. First subset consists of the vector - 
functions "rapidly rapidly" decreasing to zero for | x \ — > oo. This subset is defined by the inequality 

(6.5) \\u\\ L2 < fl-IMIc 
Here B is a fixed constant. 

Second subset consists of the vector - functions "rapidly" decreasing to zero for | x | —> oo. This subset 
is defined by the inequality 

(6.6) \\u\\ c * < B'-\\u\\ L2 
Here B' is a fixed constant. 

NIc 2 = Eilill u i]|c= = ^=1^1=0 max x | Di Ui (x) | (N = 3). 
For first subset by using appropriate sequence of inequalities [22], (|6.4p . (|6.5|) 



(6.7) ||u|| i2 < B-\\u\\ C 2 < K-\\u\\ L2 < K-\\u°\\ L2 < K-B-\\u°\ 



C 2 
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we receive 

(6.8) ||u|| c2 < K.\\vP\\ C 2, 
where K - constant. 

For second subset by using appropriate sequence of inequalities [35] , (|6.4I) , (16.61) 

(6.9) \\u\\ C 2 < B'-\\u\\ L2 < B'-\\vP\\ L2 
we receive 

(6.10) \\u\\ c * < B'.\\u°\\ L2 , 
Let us choose the constant V as: 

(6.11) max {K ■ \\vP\\ C 2, B 1 ■ \\u°\\l 2 } < V 
Then we have that 



(6.12) \\u v \\ C 2 < 1 

We can use the fixed point principle to prove existence and uniqueness of the solution of integral equation 
(& 

-= v 

For this purpose wc will operate with the following properties of matrix integral operator S v : 

1. Matrix integral operator S v depends continuously on its parameter t £ [0,oo) (based on formulas 

(15301) - @32)). 

2. Matrix integral operator Sy maps vector- functions u v from perfect space onto perfect space 
TS. This property directly follows from the properties of Fourier transform |22j , and the form of integrands 
of integral operators Sy, B (based on formulas (|3.30[) - (|3.32I0 . Hence in this case we can consider the 
convergence of functions from not only in countable system of norms ()4.2j) . but also in norm C 2 . 



3. Matrix integral operator S v is "quadratic". 

- v - V / / 

4. \\S V -U v - Sy -U V || C 2 < \\U V ~ Uy \\c 

(based on properties 1, 2, 3 and formulas p.30|) - p.32p . (|6.12p ). 



4. \S V ■ u v — S v ■ u v || C 2 < ||u v — w v ||c 2 f° r an y Wv, ^ v e TS (u v ^ u v ) and any t € [0,oo) 



-E- v 

Properties mentioned above allow us to prove that matrix integral operator Sy is a contraction operator. 



Theorem 1. Contraction operator. [23J 

Matrix integral operator § v maps perfect space T$ onto perfect space and for any u v , u v G TS 

(uy 7^ u v ) the condition 4 is valid. 

^ v 

Then matrix integral operator S v is a contraction operator, i.e. the following condition is true: 
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- V - V ) / 

(6.13) \\Sy -Uy - Sy ■ Uy\\c 2 < & • \ \ U y - U v \\c 2 

where a < 1 and is independent from u v , u v S TS for any t e [0,oo). 
Proof by contradiction. 

Let us assume that the opposite is true. Then there exist such u vn , u vn e TS (n=l,2,. . .) and 



^yn , i* VI1 
n— >oo n—too 



lim u vn , lim w,„, e TTj 

that 



.14) |[5 V • w vn - Sy ■ u vn \\ C 2 = a„ ■ ||u vn - « v „||c 2 (" = 1,2,...; a„ -» 1) 
Then the limiting result in (|6.14[) would lead to equality 

- V - V / / 

\\Sy -Uy - Sy -U V \\ C 2 = \\Uy - U V \\ C 2 , 

V . 

which is contradicting condition 4. Hence, Sy is a contraction operator. 



Theorem 2. Existence and uniqueness of solution. 

Let us consider a contraction operator S v . Then there exists the unique solution m* of equation (|6.2|) 
in space TS for any t £ [0,oo). Also in this case it is possible to obtain w* as a limit of sequence {u V n} , 
where 

- v 

Uy,n+l = S v (Uyn) (n = 0, 1, ...), 

and u v o = 0. 

The rate of conversion of the sequence {u vn } to the solution can be defined from the following inequality: 

a n 

3.15) 1 1 u vn - w v ||c 2 < 7| d|«vi ~ Uy Q \\ C 2 (n = 0, 1, ...) 

(1 - a) 



Proof: 

It is clear that 



^v,n+l = $v (^vn), Uyn = Sy [Uy^ n — l). 



It follows from (|6.13[) that 

||«v,n+l - "vnllc 2 < tt'll^vn ~ Sv,n-l||c 2 - 

Using similar inequalities one after another while decreasing n we will receive: 

||u v ,n+l - Uvn||c 2 < Ct n ■ \\Uyl - U v0 ||c 2 - 
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From this result it follows that 



\\u v , n +l ~ UvriHc 2 < ||«V,n+J ~ «v,n+i- 1 II C 2 H 1" ||#V,n+l - U vn \\ C 2 < 

(6.16) < +•••+ a n )\\u vl - u v0 \\ C 2 < - -||u v i - u v0 ||c 2 - 

(1-a) 

Because of a n — > for n — > oo the obtained estimation f|6 . 16[) shows that sequence {u vn } is a Cauchy 
sequence. Since the space T& is a perfect space, this sequence converges to an element w* G T^, 
such that S v (u* ) has sense. We use inequality (j6.13[) again and have: 

- v - v - v 

||«v,n+i - S v (u*)|| C 2 = \\Sy (u vn ) - Sy (-Sv)llc 2 < oi-\\u Vn - w*llc 2 (n = 0,1,2,...) 

The right part of the above inequality tends to for n — > oo and it means that u Vi7l+ i — > S v (u v ) 
and u* = S v (u*). In other words, w* is the solution of equation (|6.2I) . 

Uniqueness of the solution also follows from (|6 . 1 3[) . In fact, if there would exist another solution u v £ T§, 
then 

- v — - v 

\\u v - "vile 2 = ll^v (%) - S v (uv)llc 2 < a-\\uv - u*\\ C 2- 
Such situation could happen only if || u v — u^llc 2 = 0, or u v = u*. 

We can also receive an estimation (|6.15[) from estimation (|6.16[) as a limiting result for / — >• oo. 

Now let us show that continuous dependence of operator S v on t leads to continuous dependence of 
the solution of the problem on t. 

^ V 

We will say that matrix integral operator S v is continuous in t at a point to G [0,oo), if for any sequence 
{t n } € [0,oo) with t n — > to for n — > oo, the following is true: 

(6.17) S vtn (u v ) — > S vto (u v ) for any u v &TS. 

From Theorem 2 it follows that for any t G [0,oo) equation (|6.2j) has the unique solution, which depends 
on t . Let us denote it as u* t . We will say that solution of equation (|6.2p depends continuously on t at 
t = to, if for any sequence {t n } G [0,oo) with t n — > to for n oo, the following is true: 

u vt„ — * M vt • 

Theorem 3. Continuous dependence of solution on t. |23j 

Let us consider operator S vt that satisfies condition (|6.13[) for any t G [0,oo), where a is independent 

- v 

from t and that operator S vt is continuous in t at a point to G [0,oo). Then for t = to the solution of 

equation (|6.2[) depends continuously on t. 

Proof: 

Let us consider any t G [0,oo). We will construct the solution u vt of equation (|6.2j) as a limit of sequence 
{u vn } ■ 

- v 

(6.18) u v ,n+i = S vt (u vn ) (n = 0, 1, u v o = Kto) 



16 



A. TSIONSKIY, M. TSIONSKIY 



Let us rewrite inequality (|6.15[) for n = 0: 

1 

(!-«)' 



(6.19) - uvollc 2 < - — rrll^vi - u v0 \\ C 2 



Since w* to = S vto (u^ to ), then because of (|6.18j) and (|6.19|) we have: 

1 1 -- v -- v 

(6.20) ||«; t - u* v t \\c 2 < ( X _ Q ,) H" vl ~ "voiles = ^ _ ^ \\Svt [KtoJ - ^K^Hc 2 

Now with the help of (|6Tf|) we obtain the required continuity of u vt for t = to . 

Following (|6.3p and (|6.11|) we get the result: 

(6.21) u = u v • V, v = v v ■ V 2 . 

Then vector-function grad p £ T$ is defined by equations (|2.1[) where vector-function u is received from 
equation (|6.21[) . Function p is defined up to an arbitrary constant. 

In other words there exists the unique set of smooth functions Uooi(x, t), poo(x, t) (i = 1, 2, 
3) on R 3 x [0, oo) that satisfies (|27i]) . (|2~2j) . ([273]) and 

(6.22) Uooi , Poo G C°°(R 3 x [0, oo)), 
Then, using the inequality ||u||z a < ||u°||i 2 from (|6.4p . [I], [2], we have 



3.23) / |u OC; (x,t)| 2 dx<C 

JR 3 



for all t > 0. 

^ v 

Let us consider v — > in integral operator 5V ■ Then we see that theorems 1-3 arc correct also in case 
of Euler equations, i.e. there exists unique smooth solution in all time range for this case too. 

Appendix A. 

The Fourier integral can be stated in the forms: 

-| pOO pOO pOO 

U( 11 , 7 2,73) = F[u(x U X2,X 3 )} = 7 — m / / u{ Xl ,X 2 ,X 3 )e d Xl dx 2 dx 3 

(2iry/ Z J-ooJ-OoJ-OO 
-| pOO pOO nOO 

(A.l) u(x u x 2>X3 ) = 7 — m / / f/( 7 i,72,73)e -^+^-+^d 11 d l2 d l3 

(ZTT) J -oo J -oo J -oo 



The Laplace integral is usually stated in the following form: 

/•oo -i pc+ioo 

(A.2) U®{ri)=L[u{t)]= I u{t)Q-^dt u{t) = —I U®{r})Q ^drj c> c Q 

JO J c -ioo 
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(A.3) L[u (t)}=r]U®{r))-u(0) 



The CONVOLUTION THEOREM A.l. |34j . [35 
If integrals 



U®(r])= ui(t)e ~*dt Uf{n)= «a(t)e ""'di 

Jo 

absolutely converge by i?e?y > o^, then U®(r]) = Uf{rj) Ufir}) is Laplace transform of 

t 



(A.4) u{t) = / ui(t - t)u 2 {t)<1t 

Useful Laplace integral: 



o 



(A.5) 



L[e " fct ] = / e -("-».) *dt = - (Re V > r] k ) 

Jo [V ~ Vk) 



(A.6) 

In calculations of integrals (|3.30p - (|3.32[) for components of velocity u\,U2, U3 for the inverse Fourier trans- 
forms, we have each integrand /(71,72,7s) as a product of functions x(7i, 72,73) and ^(71,72,73) , 

/(7i,72,73) = x(7i,72,73) • fill, 72,7a), 
where ^(71, 72, 73) belongs to space S (functions of 71, 72, 73)[22] and x(7i, 72, 73) is one of the fractions: 

(72 +73) (71-72) (71-73) 

(7?+7 2 a +7f) ' W+7^+7!) ' (ji+lZ+ji) ' 

(73+7?) (72-73) (7 1 2 +72 2 ) 

(tT+tI+tIT' R+^Tt?)' (tT+tI+^T 

These fractions are infinitely differentiablc functions for 71 7^ 0, 72 7^ 0, 73 7^ with one point of 
discontinuity 71 =0, 72=0, 73=0. (Discontinuities have finite values at this point.) 

For these calculations the inverse Fourier transforms are defined as Lebesgue integrals with Cauchy 
principal values. 

Theorem : 



Let us prove that the inverse Fourier transform of 7(71,72,73) = x(7i, 72,7a) ■ ^(71,72,73) 
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F[f] = = / e ^^ X (l) • <Pil)dl {A.Q.I) 



ip(a) as function of a belongs to space S (functions of a). I.e., ip(cr) has two criteria: 

1) ip{ a ) is infinitely diffcrcntiable function, 

2) when a \ — > oo VK ") tends to 0, as well as its derivatives of any order, more rapidly than any 
power of jiy. 

Proof: 

Integral in (A. 6.1) admits of differentiation with respect to the parameter <jj, since the integral obtained 
after formal differentiation remains absolutely convergent: 

= r i • 7i • e ^^ X {l) ■ V(l)dl (A6.2) 

aa j J -oo 

The properties of function y>( 7 ) permit this differentiation to be continued without limit. This means 
that the function if) (a) is infinitely differentiable (see criterion 1). 

To prove criterion 2 we have created a function with parameter n: 

i 

/n(7i)72,73) = n-e ™ 2 <-'i+^+V • x(7i,72,73) ■ <p(7i7 72,73) 



n = 1,2,3... 

/n(7i>72)73) belongs to space S (functions of 71,72,73)- 

Then we estimate the integral in formula (A. 6.1) using the inverse Fourier transform of f n {l) for | a \ » 0, 
I a I — > 00 : 

/oo />oo 1 

e . x ( 7 ) . ^ l)dl \<\ n - e ^ ■ e ~~ ■ x ( 7 ) • ^( 7 )d 7 I (A6.3) 
-00 J —oo 

n » 1... < 00 

Since /n (71 7 72, 73) belongs to space S (functions of 71,72,73) then the inverse Fourier transform of 

fn (71 , 72 , 73 ) belongs to space S (functions of 01, 02, 03). Hence F[f n ] is the function, such that when 

I cr I — > 00 this function tends to 0, as well as its derivatives of any order, more rapidly than any power of 
1 

So we got from formula (A. 6. 3) that ip{cr) is the function, such that when a \ — > 00 this function 
tends to 0, more rapidly than any power of r^r. 

Then we estimate the integral in formula (A. 6. 2). Using formula like (A. 6. 3) for d Q^' , we obtain 

that 9 g^/ tends to more rapidly than any power of A. The properties of function 99(7) permit this 
differentiation to be continued without limit and further using of formulas like (A. 6. 3) to derivatives leads 
to the conclusion that all derivatives of function ip( (J ) tend to more rapidly than any power of j^. 
We have proved that V'( (7 ) belongs to space S (functions of a). 
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